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ABSTRACT 

Rotation modulates turbulence causing columnar struc- 
turing of a turbulent flow in case of sufficiently strong 
rotation. This yields significant changes in the flow charac- 
teristics and dispersion properties, which makes rotational 
turbulence modulation particularly relevant in the context 
of atmospheric and oceanic flows. Here we investigate the 
canonical flow of turbulence in a periodic box, subjected 
to rotation about a fixed vertical axis. As point of refer- 
ence we consider direct numerical simulations of homoge- 
neous isotropic turbulence. Modulation due to rotation at 
various rotation rates (i.e., different Rossby numbers) is in- 
vestigated. Special attention is paid to the alteration of 
intermittency, which is measured in terms of changes in the 
scaling of the structure functions. A reduction of intermit- 
tency quantified with the longitudinal structure functions in 
the direction perpendicular to the rotation axes will be pre- 
sented. These numerical findings correspond well to recent 
results obtained in experiments by Seiwert et al. (2008) [l]. 

TURBULENCE AND ROTATION 

Turbulence exhibits intense bursts of vorticity and strain 
that can be important for example in production of forceful 
vortices in atmospheric fiows. An important interest in at- 
mospheric fluid dynamics is concentrated around the impact 
of the Coriolis force on turbulence that tends to the two- 
dimensionalization of the flow. In this paper we consider 
direct numerical simulations of decaying turbulence in a ro- 
tating frame of reference studying the effect of the Coriolis 
force on turbulence. For high rotation rates the Coriolis force 
is dominant in a wide range of scales and plays an impor- 
tant role balancing the convective nonlinearity and viscous 
forces [2] . We will show that the Coriolis force not only sup- 
presses the forward energy transfer to small scales, but also 
modifies the dynamics of turbulence measured in terms of 
the structure functions. These are explored in this paper via 
direct numerical simulations at various rotation rates. 

The Kolmogorov K41 description of turbulence [S] |4] 
results in scaling laws for structure functions of the ve- 



locity increments. The second order structure function is 
the best-known, characterized by the famous Kolmogorov 
energy spectrum with a —5/3 slope. The Kolmogorov ap- 
proach predicts a linear dependence of the scaling exponents 
on the order of the structure function [3]. However, in 
three-dimensional isotropic turbulence a so-called anoma- 
lous scaling of the structure functions is observed [U O |6] . 
This is visible in a nonlinear dependence of the scaling expo- 
nents on the order of the structure function. This anomaly is 
associated with the effect of 'intermittency'. One may expect 
that rotation, which induces a 'trend' toward partial two- 
dimensionalization of the flow, will reduce intermittency and 
thereby also the anomalous scaling. Recent experiments in a 
freely decaying rotating turbulence using PIV show a strong 
increase of the exponents of the structure functions [l] . This 
is particularly pronounced for the second-order structure 
function. Correspondingly, a reduced scaling anomaly was 
reported. The main aim of this work is to complement these 
experimental findings with numerical simulations, allowing 
a direct correlation between reduced scaling anomaly and 
rotational flow structures. 

The organization of this paper is as follows. First, we 
introduce the computational setting for simulations of tur- 
bulence in a rotating frame of reference. Then, we present 
results of direct numerical simulations, quantifying the sup- 
pression of the energy decay for growing rotation rates. Af- 
terwards, we quantify the intermittency via computation of 
the structure functions. The tendency to two-dimensionalize 
due to rotation is clearly expressed in a reduced presence of 
smaller scales in the flow and correspondingly increases scal- 
ing exponents. Finally, concluding remarks are collected in 
the last section. 

COMPUTATIONAL SETTING 

The decay of turbulence with an additional Coriolis force 
is investigated in a simple temporal setting using a paral- 
lelized, fully de-aliased pseudo-spectral method to simulate 
the flow in a computational box endowed with periodic 
boundary conditions. The incompressible Navier-Stokes 



equations for velocity v(x, t), pressure p(x,t), and external 
force f: 



(9tv - j/Av + (v ■ V)v + Vp = f 
V ■ V = 



(1) 



are transformed into equations for the individual velocity 
Fourier coefficients u(k, t) in the wave-vector space k: 



(dt + uk^^ 



u - DW 



(2) 



In the above equation, pressure is removed using the incom- 
pressibility condition via the projection operator -Dq^ = 
<5a/3 ~ kakp/k'^; the nonlinear term W = — jr^(v • V)v) 
is transformed to the Fourier space, but solved in 

the pseudo-spectral representation; rotation is included via 
the external force also represented in the Fourier space as 
F = ^(f). Details of the code that is used for the numerical 
simulations can be found in [?]• 

Rotation involves inclusion of the Coriolis and the cen- 
trifugal force. Assuming constant rotation rates, the cen- 
trifugal force can be incorporated in the pressure term and 
the Coriolis force reads f = — 2f2 X v. In our simulations 
the rotation vector f2 is assumed to be directed along the 
z-axis; f2 = [0,0,0]. Using the Levi-Civita permutation 
symbol e and the Fourier representation, the external force 
for the rotation term can be finally written as: 



-2D„F,ne 



(3) 



In the numerical setup, the Navier-Stokes equations (|2j in 
the presence of constant rotation rate Q along the z-axis are 
reformulated using the complex helical-wave decomposition 
framework '8': 



dt + vk^ I 



fc2 



Q^AQ 



Q-1u = Q-1dW, (4) 



where I is the unit matrix. The matrix A is defined as: 
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and the matrices Q and Q~^, corresponding to the eigen- 
values of the matrix A: Ai = +ik, A2 = —ik, A3 = after 
solving the characteristic equation, have the form: 
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(7) 

where q = (2fe^s)~^/-^ and s = + The Coriolis force 
is diagonalized through this transformation and introduced 
as an integrating factor along with the viscous term in the 
Navier-Stokes equations. For time-integration the fourth- 
order, four-stage, Runge-Kutta method is used. To remove 
the aliasing errors, the nonlinear term is computed on two 
grids (original and shifted) with an additional truncation 
to m = 120 active modes for the lower resolution (256'^) 
and m = 241 for the higher resolution (512^) runs at the 
cut-off number kc = 27r(m -|- 1/2) /Li, for the cubic box of 
length L5 = 1. Diagonalization of the Coriolis force through 



helical wave decomposition allows performing simulations 
at acceptable time-steps as the rotation does not dominate 
the reduction of the step in the explicit integration of the 
equations l|4}. Flow results obtained with direct numerical 
simulation are described next. 




Figure 1: Vertical structuring of rotating turbulence at Rx = 
92. Snapshot of vorticity in j/-direction at t = 2.5 for Ro = 
00 (a) and Ro = 0.05 (b). 



DIRECT NUMERICAL SIMULATIONS 

Direct numerical simulations of the spin-up of ini- 
tially homogeneous isotropic turbulence at various rotation 
rates Q and an initial Taylor-Reynolds numbers Rx = 
50, 100, 200 were performed. The decay of turbulence in 
the rotating setting at Rossby numbers Ro = = 
00,0.5,0.25,0.1,0.05,0.025,0.01,0.005, which correspond to 
O = 0, 1, 2, 5, 10, 20, 50, 100, 200 at reference velocity Ur = 1, 



is systematically investigated. The initial condition for these 
simulations was taken from a simulation of non-rotating, ho- 
mogeneous, isotropic decaying turbulence at a resolution 
of 512''. The generation of this condition is based on a 
random velocity field with prescribed spectrum, which is 
decayed for approximately two initial eddy-turnover times 
TO- (t = 0.5 ~ 2to, where tq denotes the initial eddy- 
turnover time). Subsequently, the evolved field is rescaled 
to the initial energy value (Eq = 0.5n^) to serve as the ini- 
tial condition in the simulations for various rotation rates. 
This procedure preserves the initial Reynolds number. Sim- 
ulations at various Rossby numbers were performed with a 
resolution of 256^ for the lower Reynolds number and 512'^ 
for higher turbulence levels. We simulated up to t = 6 
(« 24ro). 

A snapshot showing the component of vorticity in the 
i/-direction without rotation is presented in Fig. [Ha). The 
small-scale vortical flow structures are seen to undergo strik- 
ing qualitative changes when the system is rotating about 
the ^-axis, see Fig. [TJb) representing the flow at Ro = 0.05. 
The rotational flow structuring corresponds to a trend to- 
ward (quasi) two-dimensional flow reminiscent of a strict 
Taylor-Proudman state. 
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Figure 2: Decay of the energy E{t) for various rotation rates 
(Rossby numbers) at the initial R\ = 200. Curves are la- 
beled with Ro. 

Corresponding to the much larger coherent structures in 
the flow, the decay rate of the kinetic energy E is strongly 
reduced due to rotation as may be observed in Fig. [2] The 
energy decay is usually characterized by a decay exponent n 
such that E{t) ~ t" for large times [9l llOj. In case of self- 
similar decay consistent with E{t) ~ t", the decay exponent 
n should be constant in time. Hence, the results presented 
in Fig. [2] should form straight lines when shown on a log- 
log scale. The decay exponent associated with the slope of 
the decay can then be easily examined. In our simulations, 
two separate regions can be associated with the decay of the 
energy. Only for sufficiently large times (see Fig. [Sj do we 
observe self-similarity. In the first period up to about t = 1, 
the energy decay slope significantly varies in the first period 
of the decay. This is in particular visible for low rotation 
rates. Such behavior can be understood from the fact that 



at low rotation rates the flow needs more time to adapt itself 
to the imposed rotation. The isotropic homogeneous turbu- 
lence is preserved for a longer time in such cases. For high 
rotation rates, all flow scales are strongly influenced and a 
transitional period before reaching the self-similar stages is 
much less pronounced. Consequently, the self-similar decay 
of energy can be safely examined for t > 1. The structure 
functions were also computed for t > I. The decay exponent 
is observed to decrease as the rotation rate is increased. Such 
behavior of the decay exponent is in agreement with results 
obtained for the experimentally investigated grid-generated 
turbulence in a rotating tank During the initial stage 

of the decay, the energy decay is more suppressed for higher 
rotation rates and consequently the energy dissipation rate 
has much smaller values in the initial phase of the decay, 
see Fig. |4] Similarly, as in case of the kinetic energy, the 
self-similar decay of the energy dissipation rate is present 
for large t only. This is particularly true for the smaller 
rotation rates. 
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Figure 3: Log-log plot of the energy decay presented in 
Fig. El 

To assess the resolution of the flow we examine the so- 
called femaxf? criterion, where rj is the Kolmogorov length 
and fcmax is the largest wavenumber used in the simula- 
tions. The result is shown in Fig. |5] In order to resolve 
dynamically relevant turbulence length-scales, it is required 
that femax»? is sufficiently large. A commonly accepted cri- 
terion of adequate spatial resolution is that fcmaxf? > 1- For 
all rotation rates, this criterion is satisfled. This underpins 
confidence that the Kolmogorov scale is properly resolved 
in our simulations. Larger anisotropic structures built by 
backward energy transfer for larger rotation rates. As the 
flow develops in time, the energy decays further and the 
energy dissipation rate decreases to very small values. Sub- 
sequently, quite large values of the Kolmogorov scale rj can 
be observed for the final periods of the decay. 

Rotation changes the statistical properties of the turbu- 
lent flow. For example skewness is reduced from its theoret- 
ical value oi S = 0.5 in homogeneous isotropic turbulence as 
proven by Batchelor |12| . cf. Fig. |6] In our simulations we 
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Figure 5: Resolution kma-KV criterion for various Rossby 
numbers as a function of time at the initial = 200. 

use the spectral formula for the computation of the skewness: 



■^^'^ 35 («(*)) 



k^T{k,t)dk, 



(8) 



where X(t) is the Taylor microscale, u{t) is the rms ve- 
locity, and T(k,t) is the energy transfer term at time t 
and wavenumber k. The observed significant modulation 
of skewness for higher rotation rates can also be associ- 
ated with a departure from isotropy in the considered flow. 
The current simulations display an immediate acceleration 
in the initial period of the flow development. This imme- 
diate change in the flow properties for high rotation rates 
is seen through a significant reduction of skewness in the 
first time-steps of the simulations. These initial data are 



Figure 6: Skewness S for various rotation rates (Rossby 
numbers) as a function of time at the initial Rx = 200. 



INTERMITTENCY 

The structure functions used in our setup have the fol- 
lowing form: 

Sj,{r) = (\5vrn , (9) 

where p is the order of the structure function and &Vr = 
(v(x+r) — v) r/r is the longitudinal velocity increment com- 
puted as the volume average (■). The separation distance r 
normal to the rotation vector f2. Intermittency usually refers 
to the anomalous scaling ?p in the similarity hypothesis ap- 
plied for the structure functions: 



Sp{r) oc r'^p, 



(10) 



The scaling exponents were computed from a numerical 
derivative: = (dlog(5p(r)))/((i log r). 

Theory of homogeneous isotropic turbulence predicts 
that in a self-similar flow an inertial range exists in which 
?p = p/3 13 and correspondingly, e.g., the energy spec- 
trum is described as E(k) ~ k~^^^ . In contrast, for rotating 
turbulence much steeper energy spectra are predicted by 
Kraichnan assuming a totally inhibited energy transfer. This 
was predicted to yield E{k) ~ k~^ or E{k) ~ k~^ for a 
spectrum dictated by an energy transfer time-scale propor- 
tional to l/C Il4l . Subsequently, the anomalous scaling of 
the structure function = p/2 was experimentally found 
in rapidly rotating fluids 1151 . Reduced intermittency mea- 
sured as a departure from the ?p/?2 = p/2 will be shown 
here for the rotating decaying turbulence. 

An example of the structure functions for increasing or- 
der p is shown in Fig. [7] (compare with Fig. 2 in 1 ). The 
separations r are normal to the rotation vector O. This is 
particularly important, as the observed departures from the 
isotropic properties of turbulence are analyzed in a direction, 
which is perpendicular to the rotation axis. Hence, they are 



not directly influenced by the anisotropic character of the 
flow in this direction. The presented investigations are re- 
stricted to the self-similar range t > 1 and compared to the 
ESS dSlin] method. 





(b) 

Figure 7: Structure functions for increasing order at an ini- 
tial R\ = 200 at t = 2 and rotation rate £7 = 1 (Ro = 0.5) 
plotted as a function of the separation scale r (a) and as 
a function of 52 (b). The curves for p = 4,6,8,10 have 
been vertically shifted by factors 10^,10**, 10^ and 10* for 
visibility. 

Our DNS results show modulation of the structure func- 
tions due to rotation, as illustrated in Fig. [7] for a particular 
instance in time and low rotation rate (Ro = 0.5). The raw 
?p and normalized exponents fp/<;2 at different times during 
the decay are shown in Fig.[8] Similarly, as in the experimen- 
tal results [l], the exponents differ from the strictly linear 
p/2 law reported in |15| . 

The reduction of intermittency can be clearly seen from 
the intermittency factors -jp = p/2— fp/<;2, which should van- 
ish for non-intermittent velocity fluctuations. These factors 
were plotted for various instances in time and low rotation 
rate in Fig. [9] They show similar behavior as in the experi- 
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Figure 8: Structure functions exponents (a) and normalized 
exponents fp/?2 at an initial Rx = 200 at various times 
during the decay used to label the curves, using = 1. 
Scaling exponents obtained from She-Leveque (a) and p/2 
(b) are indicated with the thick line. 

ment with a significant decrease during the initial stages of 
the flow development. This behavior is probably connected 
to the initial rapid spin-up of the flow that causes immedi- 
ate two-dimensionalization of the flow, i.e., the formation of 
elongated flow structures along the axis of rotation. 

SUMMARY 

Our numerical simulations of rotating decaying turbu- 
lence indicate a strong increase of the structure function 
exponents during the decay. The observed modulation of 
the structure functions is directly associated with changes 
in the flow characteristics, i.e., the increase of the structure 
function exponents is related to the much steeper energy 
spectra in case of rotation. An important fact is that the 
structure functions are measured in a direction perpendicu- 
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Figure 9: Intermittency factors 7p = p/2 — fp/<;2 for large 
Rossby number Ro = 0.5. 

lar to the axis of rotation. This implies that the observed 
changes of the turbulence characteristics appear in a di- 
rection, which is not directly influenced by the anisotropy 
introduced through the Coriolis force. Non-negligible devi- 
ations from Kolmogorov-universality and self-similarity are 
observed as a result of anisotropization of the flow. This de- 
parture can be directly linked with the modulation/interplay 
of the intermittency phenomenon and scale-dependent en- 
ergy transfer in the energy cascade process (deviation from 
homogeneity). As a result, it may be inferred that in case 
of rotating turbulence, both the longitudinal and transverse 
components of the structure functions are equally relevant, 
as the contribution to the energy transfer from the latter 
component does not statistically vanish because of homo- 
geneity as in the K41/K62 framework. In particular, the 
presented results for rotating decaying flow are in line with 
experimental findings obtained in a rotating tank with back- 
ground turbulence [I] and forced flow conditions obtained 
through direct numerical simulations 1181 . Strong increase 
of the second order scaling exponents suggests that in case of 
rotation, steeper energy spectra are predicted than those ob- 
tained under the assumption of nonlinear interactions lead- 
ing to an energy transfer governed by the time-scale 
E{k) ~ 1141 . Further research will be devoted to this 

subject, in which energy spectra exponents will be analyzed. 
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